Astronomy & Astrophysics manuscript no. 


February 2, 2008 


(DOI: will be inserted by hand later) 





Semi-classical theory of collisional depolarization of spectral lines 
by atomic hydrogen. I: application to p states of neutral atoms. 

M. Derouich 1 , S. Sahal-Brechot 1 , P. S. Barklem 2 , and B. J. 0'Mara t3 



o 
o 

(N 
c3 



in 
> 
in 



o 
m 
o 

S3 

Oh 

6 

03 



X 



1 Observatoire de Paris-Meudon, LERMA FRE CNRS 2460, 5, Place Jules Janssen, F-92195 Meudon Cedex, 
France. 

2 Department of Astronomy and Space Physics, Uppsala University, Box 515, S 751 20 Uppsala, Sweden 

3 Department of Physics, The University of Queensland, St Lucia 4072, Australia 
e-mail: Moncef . Derouich@obspm . f r 

Received 2002 / accepted 2003 

Abstract. The present paper extends the method of Anstee, Barklem and O'Mara (Anstee 1992 Anstee & 
Q'Maxa lT55nH555l Anstee, O'Mara & Ross tlWl Barklem HMfl Barklem & O'Mara H557I Barklem, O'Mara & 
Ross ll99*5l . developed during the 1990's for collisional line broadening by atomic hydrogen, to the depolarization 
of spectral lines of neutral atoms by collisions with atomic hydrogen. In the present paper, we will limit the 
calculations to p (I = 1) atomic levels. The depolarization cross sections and depolarization rates are computed. 
In Table |21 cross sections as functions of the relative velocity and effective quantum number are given, allowing 
for the computation for any p atomic level. Our results are compared to quantum chemistry calculations where 
possible. The sensitivity of depolarization cross sections to regions of the potential is examined. We conclude that 
the accuracy obtained with our method (< 20 % for the depolarization rates) is promising for its extension to 
higher i-values for the interpretation of the "second solar spectrum". This will be the object of further papers. 
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1. Introduction 

Observations of linear polarization at the limb of the Sun 
produce rich structures. This polarization spectrum (the 
so-called "second solar spectrum" ) is an entirely new spec- 
trum to be explored and interpreted (Stenfo & Keller 
1997). This polarization is due to scattering of the incident 
anisotropic solar radiation field. The radiative anisotropy 
is due to the structure of the semi-infinite or layer atmo- 
sphere in the vicinity of the surface, and is revealed by the 
center-to-limb darkening of the solar brightness. 

In order to interpret quantitatively the Stokes pa- 
rameters of the observed lines, a non-optically-thin and 
non-LTE model in the presence of a magnetic field has 
to be used. In fact, collision induced transitions between 
Zeeman sublevels can play a major role in depolarizing 
the levels (and thus the lines). Depolarizing collisions of 
the radiating and absorbing atoms with neutral hydrogen 
atoms of the medium have to be included in the statis- 
tical equilibrium equations. Very few depolarization rates 
have currently been computed, they have been obtained 
through sophisticated quantum chemistry methods which 
are accurate but cumbersome. It would be useful to de- 
velop methods capable of giving reasonable results for 
many levels rapidly. This is the object of the present work. 



During the 1990's Anstee, Barklem and O'Mara 
(Anstee Anstee & O'Mara HMTHT^l Anstee et al. 
ITW71 BarklemEHBI Barklem & O'Mara HWI Barklem et 
al. I1998f) developed a new semi-classical theory, in order 
to obtain spectral line widths due to collisions with atomic 
hydrogen that were better than those obtained by the Van 
der Waals approximation. In this theory, which will be 
referred as ABO in the following, the scattering S ma- 
trix is obtained by solving semi-classical coupled differen- 
tial equations. The hydrogen-atom interaction potential is 
derived from time-independent second-order perturbation 
theory without exchange, allowing the Lindholm-Foley av- 
erage over m states (Brueckner 119711 O'Mara I1976[) to 
be removed. Whereas the Van der Waals potential un- 
derestimates the broadening by typically a factor two or 
more, the ABO method gives rather good agreement (< 
20 % for the spectral line widths) with the results obtained 
from quantum chemistry calculations. In addition, these 
authors used their ABO method by applying it to a se- 
lection of strong solar lines. The derived abundances are 
consistent with meteoritic ones. 

In the present paper, we extend the ABO theory to 
the calculation of collisional depolarization rates which 
enter the statistical equilibrium equations for the polar- 
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ized atomic density matrix. This paper is focused on the 
method and is applied to p states of neutral atoms. We will 
present and compare our results with those recently ob- 
tained by quantum chemistry methods. We also compare 
the results to the depolarization rates calculated with the 
Van der Waals interaction. A discussion of our results is 
given, studying the region of the potential curves which 
play the major role in the variation of the depolarization 
cross sections. In conclusion, we show that our method, 
owing to the simplicity of the processes taken into account, 
is a useful alternative to the accurate but time consuming 
quantum chemistry method for obtaining most of the de- 
polarization rates to better than 20 % accuracy. It should 
now be possible to rapidly obtain the large amount of data 
needed for the interpretation of the second solar spectrum, 
in particular for the heavy atoms that are inaccessible to 
present quantum chemistry calculations. An extension to 
d and / states is in progress. An extension to ions is in 
progress and this will be the subject of further papers. 

2. Formulation of the problem 

2.1. The statistical equilibrium equations (SEE) for 
polarization studies 

The calculation of the Stokes parameters of a spectral 
line observed at the solar limb is a non-LTE (Local 
Thermodynamical Equlibrium) problem. It needs the solu- 
tion of the coupling between the polarized radiative trans- 
fer and the statistical equilibrium equations (SEE) for the 
populations and coherences of the polarized atomic den- 
sity matrix. In fact the atomic levels are only "aligned", 
which means that the observed lines are globally linearly 
polarized. This linear polarization is due to anisotropic 
scattering of the incident solar radiation field which con- 
tributes to excitation of the atomic levels, the anisotropy 
being due to limb darkening, even at zero altitude above 
the limb. 

Sahal-Brechot l|1977[) developed the SEE for an 
"aligned" multilevel atom, i.e. for the Zeeman sub- 
levels \a,JMj). She took into account excitation by the 
anisotropic solar radiation field, deexcitation by sponta- 
neous and induced emission, and excitation, deexcitation, 
polarization transfer and depolarization by isotropic colli- 
sions of the particles of the medium. She included the ef- 
fect of a "strong" magnetic field (strong case of the Hanle 
effect) and thus coherences between Zeeman sublevels 
could be ignored. Only the populations of the Zeeman 
sublevels were taken into account. 

She wrote these equations both in the dyadic basis 
\a JMj) (a' J'M'j\ built on the standard basis of the 
atomic states \aJMj) and in the irreducible tensorial 
operator basis aJa J T*. The expansion of the aJ a J 
over the dyadic basis is given by Fano & Racah 1)1959)) 
(cf. also Messiah ITMT1 Fano 11555) 1: 

aJa'J'rpk _ (_-\\j'—M'j 

x \aJMj)(a'J'M'j\(JJ'Mj - Mj\k - q), (1) 



where (JJ'Mj — M'j\k — q) is a Clebsch-Gordan coeffi- 
cient. She demonstrated the interest and power of this 
second basis, which is especially well adapted to problems 
with spherical and/or cylindrical symmetry processes. The 
equations are simplified, and the physics is extracted and 
highlighted. In particular, the odd fc-terms can be elimi- 
nated from the SEE when the incident radiation field is 
anisotropic but unpolarized. In fact the natural radiation 
field transports k = (population) and k — 2 (align- 
ment) terms only. Orientation terms (k = 1) vanish in the 
present astrophysical problem. 

Sahal-Brechot l|1977[) gave the general expressions for 
the various transition probabilities due to isotropic and 
anisotropic radiation and for isotropic collisions, i.e. ex- 
pressions for population, orientation and alignment trans- 
fer between the levels, creation and destruction of popu- 
lation, orientation and alignment (depolarization rates in 
particular). At the same time, Omont 1)1977(1 developed a 
general theory of relaxation for a two-level atom. In addi- 
tion he introduced coherences. 

In fact, for isotropic processes, the coupling terms in 
the SEE implying transfer of coherence q — > q' are zero, 
and due to the isotropy of the collisions with hydrogen, 
the depolarization rates are g-independent. This is one of 
the advantages of the irreducible tensorial operator ba- 
sis. Therefore, as we will consider isotropic collisions only 
in the following, it will be sufficient to calculate colli- 
sional cross sections between the sublevels, as for the no- 
coherence case ( Sahal-Brechot IT977[I . even if coherences 
play a role in the SEE. We recall that coherences can ap- 
pear when the magnetic field is weak and not directed 
along the solar radius (Hanle effect). We also notice that 
in the present paper we are interested in hydrogen-atom 
collisions and not in electron or proton-ion collisions as in 
Sahal-Brechot 1)1977)1 . After integration of the cross sec- 
tions over a Maxwell distribution of velocities and multi- 
plying by the local hydrogen density, the depolarization 
rates obtained (and eventually polarization transfer rates 
between levels) can be used in the statistical equilibrium 
equations. 

We focus on the collisional depolarization rates in the 
present paper. We will use both the dyadic and the irre- 
ducible tensorial operator bases. The latter one is in fact 
more convenient. 

An important preliminary point to notice is the fact 
that the radiating atom can have a nuclear spin and thus 
hyperfine structure. The importance of the effect of the 
hyperfine structure on the Stokes parameters of the line 
studied was discussed for the example of the hydrogen 
coronal line Lya by Bommier & Sahal-Brechot 0982). 
They showed that the SEE must be solved for the hy- 
perfine levels when the inverse of the lifetime of the upper 
level is smaller than the hyperfine splitting, i.e. the hyper- 
fine levels are separated. This arises from the Heisenberg 
uncertainty principle. On the contrary, when the hyperfine 
splitting is small compared to the inverse of the lifetime, 
the hyperfine levels are mixed and degenerate, the hyper- 
fine structure can be ignored and the SEE can be solved 
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F 


AEhfs('P 3 /2) 


AE HF s( 2 S 1/2 ) 


3 


58.44 




2 


34.32 


1772 


1 


15.86 





Table 1. Values of the energy separation between the hy- 
perhne levels F and F — 1 for: Na in the ground state, 
3s 2 Si/ 2 (Ackermann I1966[) . and in the excited state 
3p 2 P 3 / 2 (Gangrsky et al ll998[l . All energies are in MHz. 



for the | a J) fine structure levels. This was the case for the 
hydrogen coronal line where the F levels are degenerate 
and the J levels are well separated. 

In the second solar spectrum, where metallic lines are 
observed, the inverse of the lifetime (9.7 MHz for the Na 
D lines) (Gaupp et al 1982) is often small compared to the 
hyperfine splitting (if it exists) (cf. TableUJ). Then the SEE 
must be solved for the \aJF) hyperfine levels (Bommier 
& MolodijEEEl Kerkeni & Bommier EEEJ) and collisional 
rates between the hyperfine sublevels must be calculated. 

2.2. Dyadic basis for the depolarization rates 

The depolarization arises from isotropic collisions between 
the Zeeman sublevels of a level \aJ): the collisions reestab- 
lish thermodynamical equlibrium between the sublevels 
and thus atomic polarization vanishes. We have to calcu- 
late (Sahal-Brechot H977f) for each sublevel \aJMj) (or 
\aJFMp) if there is an hyperfine structure): 



D(aJMj,T) 



C(aJMj 



aJM'j,T) 



(2) 



where ^{aJMj — > aJM'j,T) is the collisional transition 
rate between the sublevels \aJMj) — ► \aJMj). It can be 
written as a function of the local temperature T and the 
hydrogen perturber local density Uh- 

((aJMj -> aJM'j, T) = n H I a(aJMj -> aJM'j, v) 

Jo 

xvf(v)dv, (3) 

f(v) being the distribution function for the atom- 
perturber relative velocity v, and a(aJMj — > aJM'j,v) 
is the cross section. For the Maxwellian distribution, 



(4) 



where p, is the reduced mass of the system. 

The relaxation term for the density matrix (diagonal 
term of the system of SEE) due to collisions is given by: 



d a PMjMj 
dt 



)coit = -D{aJMj,T) aJ p M jMj (5) 
- E ^2((aJMj^aJ'M'j,T) 

a'J'^aJ M'j 
x PMjM'j 



The terms added to the depolarization term D(aJMj, T) 
are the so-called quenching terms. 

The dyadic basis is not suitable because the depo- 
larization terms depend on the chosen dyadic basis. The 
aj,a j j,k k as } s i s we ll adapted because the depolarization 
terms do not depend on the basis and have a physical 
interpretation. 



2.3. Basis of irreducible tensorial operators 

In this basis, the SEE relaxation term in the non- 
coherence case (q = 0) fSahal-Brechot !1977l OmontHHUJl 

becomes 

^-^)coii = -D k (aJ, T) aJaJ pQ + quenching term, (6) 

where aJ Pq is the multipole fc-moment of the density ma- 
trix for the level \aJ). For isotropic collisions D k (aJ,T) 
is characterized by (2J + 1) constants (0 < k < 2J) 
and D°(aJ,T) is called the destruction rate of popula- 
tion, which is zero for the no-quenching approximation. 
D l (aJ,T) is the destruction rate of orientation (circular 
atomic polarization). D 2 {aJ,T) is the destruction rate of 
alignment which is of interest in astrophysics (linear po- 
larization), and the depolarization rate D k (aJ, T) is given 
by: 

poo 

D k (aJ,T)=n H a k (aJ,v)vf(v)dv (7) 

Using our present notation, formula (20) in Sahal-Brcchot 
H1977JI with k = k' can be written as: 



a k (aJ,v) = (2fc + l)£ J2 (- 1 ) 2 ' 7 



-2M, 



M'j Mj^Mj 

J k J 
-Mj Mj 



a{aJMj -> aJM'j,v) 



-(2k+l)J2 (-l) 2J - Mj - M ' J (8) 



Mj M'j^Mj 

J k J 



-Mj Mj ] \ -M' 7 M'j 



J k J 

J M'j 



xa(aJMj -> aJM'j,v) 

The expressions between parentheses denote 3j- 
coefficients (Messiah I1961|l . After some calculations 
which are not detailed here, this equation can also be 
written as: 

a k (aJ,v) = a(aJ,v)-(2k + l) ^ f^fJ-Mj-Wj 

M.,,M'j 

J k J \ ( J k J 
-Mj Mj ) 

xa(aJMj -> aJM'j,v) 



M'j M'j 



(9) 



Where a(aJ, v) is the elastic cross section for the level 
| a J) and the relative velocity v: 



a(aJ, v) = a(aJMj — > aJM'j, v) 



(10) 



M' 
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This cross section can also be written as an average over 
the initial states Mj as usual: 

a{aJ ' v) = 2JTT ^ cr(aJMj ^ aJM'j,v) (11) 

Mj,M'j 

where the angular average over all the directions of the 
collisions is not necessary here because the averaging over 
Mj is equivalent to the angular average over the directions 
of the collision. This property will be used in Sect. [SJ first 
line in equation (|41|l . 

Since the depolarization cross section a k (aJ,v) can 
be written as a linear combination of the a(aJMj — > 
aJM'j,v), the depolarization rate D k (aJ,v) is a linear 
combination of the collisional rates between the Zeeman 
sublevels ((a,JM -> aJM') fSahal-Brechot ll977[> . In par- 
ticular we have (Sahal-Brechot 1977), for J = 1: 

D x {al) = C(10 -> 1 ± 1) + 2C(1 =F 1 — > 1 ± 1) 
D 2 {al) = 3C(10^ 1±1), 



where I is the unit matrix and T = I — S is the so-called 
transition matrix. 

3.1. Calculation of the S-matrix 

The Hamiltonian of the system atom-perturber can be 
written as: 



H = Ha + Hp + V, 



(14) 



and for J = 3/2 (Gordeyev et al HMVl Roueff & Suzor 
1974): 



*(«§) 



5 SV 2 2 
8 3-3 

5 SV 2 2 



D 2 {a-) = 2C( 



D\a-) 



3 3 

r, >/3 -3 
2C( 2- 
9 3-1 

5 SV 2 2 
2 3-3 

5 SV 2 2 



2 2' 

2 2' 
3 1, 

2 2 ) + 
2 2 7 
2 2' 
2 2' 



2 33 
5 s ^2 2 



9 

5^2 ~2 



--) 

2 2 ; 

3-3 3 3 



2 2 



2C( 



3 -1 
2 "2" 



8 33 
5 sv 2 2 



3 1, 
2 2 ; 



3 1 



where Ha is the atomic Hamiltonian, Hp that of the per- 
turber and V is the interaction energy. Ha can be written 
as: 

H a =H +H fs + H hfs , (15) 

where H$ contains the electrostatic part, Hps the fine 
structure and Hhfs the hyperfine structure if it exists. 

It can be shown that the hyperfine splitting AEhfs 
is always small when compared to V and thus J and I 
(the nuclear spin) can be decoupled during the collision 
process. This arises from the Heisenberg uncertainty prin- 
ciple: if r is a typical collision time, one can ignore the 
hyperfine structure if 

r AEhfs « 1, (16) 

where r ~ — , b being a typical impact-parameter (of the 

v 

order of 10ao, ao being the Bohr radius) and v ~ 4.8 x 
10~ 3 a.u. for a temperature of about 5000 K so that t ~ 
2.08 x 10 3 a.u.. In the case of the Na atom the nuclear spin 
is / = 3/2, which yields two hyperfine levels in the ground 

3/2 



state 3s z S 1/2 , F = 1 and F = 2. For the excited 3p 2 P : 



2 2' 

1 3 :\ :-! :-! 



~2 2 



2 2' 



3. Calculation of the cross sections 



We will use a semi-classical treatment for the collision: 
we will assume that the hydrogen perturbing atom moves 
along a straight path, unperturbed by the collision process 
and characterized by an impact-parameter b. 

In the dyadic basis \aJMj){a' J'M'j\ (or 
\aJFMp)(a' J' F 1 M' F \ if there is an hyperfine struc- 
ture), the cross section follows from an integration over 
the impact-parameter b and averaging the transition 
probability P(aJMj — > aJM'j,b,v) over all the collision 
directions: 

a(aJMj -» aJM'j,v) = 27r6d6 f ^ (12) ( a JFM F \S\a'J'F'M' F ) = 



j> P(aJMj -> aJM'j, b, v)dQ 



and then the collisional transition rate is given by equation 
J3j. The transition probability can be written in terms of 
the 5-matrix elements, 



P(aJMj -> aJM'j, b, v) 



\{aJMj\I - S{b,v)\aJM t 



(13) 

' \ |2 



state there are four hyperfine levels: F = 0, F = 1, F = 2 
and F = 3. The energy differences between these levels are 
given in Table fj] In the case of the ground state 3s 2 6*1/2, 
AE HFS ~ 2.694 x 10~ 7 a.u. (1772 MHz) (Ackermann 
rTM>|) and then r AE HFS ~ 5.6 x 10~ 4 < 1. For the 
excited level the most important hyperfine splitting is, 
of course, between the levels F = and F = 3 where 
AEhfs - 1-653 x 10 _8 a.u. (108.72 MHz) (Gangrsky et 
allESHI) and then t AE hfs ~ 3.4 x 10~ 5 < 1. 

The hyperfine splitting can always be ignored in the 
treatment of the collision under these conditions. We have 
F = I + J and M F = Mj + Mj, and 

\aJFM F ) = (/JMiMj|FM f )|JMj)|7M/), (17) 

M u Mj 

where (IJMjMj\FMp) is a Clebsch-Gordan coefficient. 
Then the matrix element 

(IJMiMj\FM f ) 

Mi,Mj,M' t ,Mj 

x (IJ'MjM'jIF'M'f) (aJMj \ (IMj \S\IM'j)\a' J'M'j) . (18) 

If the hyperfine structure can be neglected during the col- 
lision, the S-matrix is diagonal in I and its elements do 
not depend on M/ {I is conserved). Thus M[ = Mj and 

(aJFMF\S\a'J'F'M F ) = (IJMjMj | FMp) 

Mi,Mj,M'j 

x {IJ'MiMj\F'M' f ) (aJMj\S\a' J'M'j) (19) 
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This relation holds for each impact-parameter and each 
relative velocity. 

We recall the relation between Clebsch-Gordan coeffi- 
cients and 3j-coefficients: 

(jij2m im2 \jm) = ^2j+T{-l)^-^ +m ( 31 32 3 ) 
v \ mi m,2 —m J 

We use the convention of Messiah H19fil|l for the phase 

factors. 

3.2. Calculation of the S-matrix elements between fine 
sublevels 

The fine structure splitting can generally be neglected dur- 
ing the collision problem at the temperature of interest 
in the second solar spectrum (~ 5000^). The validity 
of the condition is the same as above, except that the 
fine structure splitting is larger than the hyperfine split- 
ting. For example in the case of the Na atom, at 5000.K", 
AE FS ~ 0.78 x lCT 4 a.u. and r ~ 2.08 x 10 3 a.u and so 
t AEps — 0.16 < 1. This has been discussed in earlier 
papers focussed on collisional line broadening, see for ex- 
ample Roueff ljl974ll and related papers. She showed that 
the fine structure splitting is of importance at low temper- 
atures (200-500 K) for the sodium D lines and negligible 
at higher temperatures for line broadening studies (500 K 
and higher). The D\ and D 2 lines have in fact the same 
collisional width for hydrogen perturbers at solar temper- 
atures. 

If the fine structure is neglected, where J = S + 1, and 
Mj = Ms + mi, we have: 

(edJMj\S\a'l'J'M'j) = ^ (SlM s mi\JMj) (20) 
x{Sl l M s m l l \j'M l J ){almi\S\a'l'm' l ), 

cf. equation ljT§|) . 

However, we notice that the parallelism between line 
broadening and depolarization cross sections must be used 
carefully in the peculiar following case. Depolarization of 
hyperfine components of spherically symmetric J = J' = 
1/2 levels (ground state of sodium for instance), needs the 
matrix elements Mj = ± I — > Mj = ± A . If the fine struc- 
ture is neglected, they are zero and the depolarization and 
polarization transfer rates of the hyperfine components of 
the ground states are zero. Kerkeni et al (200C§ calcu- 
lated the hyperfine collisional depolarization rates for the 
ground state of sodium with a quantum chemistry method 
which takes into account the fine structure splitting dur- 
ing the interaction processes. They showed that these rates 
are of the same order of magnitude as those of upper levels 
for which the fine structure splitting plays a minor role in 
the collisional processes. In fact this remark does not con- 
cern I > 1 states and is very specific. However, we consider 
that it is important to mention it. 



4. The model of Anstee, Barklem and O'Mara 
(ABO) 

The present semi-classical model used for the depolariza- 
tion is the same one as that already developed for line 
broadening theory applications in the 1990's by Anstee, 
Barklem and O'Mara (Anstee Anstee & O'Mara 

ITMT1 [15551 Ans tee et al. IT5571 Bark lem IT5581 Barklem 
& O'Mara lTMTl Barklem et al. [T55g|l . 




Fig. 1. The perturbed atom core is located at A and the 
hydrogen perturbing core (a proton) at P. Their valence 
electrons are denoted by 1 and 2 respectively. 



The theory is described in these papers, and thus only 
a brief review is presented in the following: 

1. The model is of semi-classical nature, the internal 
states of the two atoms are treated within the framework 
of quantum mechanics and the collision is modelled with 
the straight path classical approximation: it is assumed 
that the interparticle interaction is too weak to deflect 
the perturbed particle from the straight classical path. 

2. The interaction energy between the two atoms is 
smaller than the energy eigenvalues of the isolated atoms, 
so that the Rayleigh-Schrodinger perturbation theory can 
be used for obtaining perturbed eigenvalues. Ionic and ex- 
change interactions are neglected. The hyperfine and fine 
structure interactions are also neglected. 

3. Time-independent perturbation theory is applied to 
second order. It is assumed that the energy denomina- 
tor in the second order term can be replaced by a suit- 
able average energy E p . This is the Unsold approxima- 
tion (Unsold E23 Unsold ES5J. The value E p = -4/9, 
in atomic units, is adopted. This approximation works 
when the energy level separations of the perturbed atom 
are small compared with those of hydrogen. For neutral 
atoms this condition is generally satisfied. This is particu- 
larly valid when considering interactions with hydrogen in 
the ground state, like in this work, which is well separated 
from any other levels (Anstee & O'Mara 1991). 

4. Unperturbed atomic eigenfunctions are used. The 
Coulomb approximation with quantum defect is used. 
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5. Treatment of the collision 

Each collision is characterized by an impact-parameter 
vector b and a relative velocity vector v. Since the hyper- 
fine and fine structure interactions are neglected, equation 
(|14|) becomes: 



H = Hn + H f 



V 



(21) 



where V is the interaction energy which is assumed to 
be totally electrostatic. Only V depends on time. In the 
coordinate system of Fig. ^ V is given, in atomic units, 
by (Anstee & 0'MaraHj991): 




Fig. 2. Trajectory of the perturber (H) in the two ref- 
erence frames. R is the interatomic separation, R = 
yjb 2 + v 2 t 2 . Axyz and AXYZ are the atomic and rotating 
frames. 



V = 



R 



1 

ri2 



1 

T2 



1 

Pi 



(22) 



and atomic units are used hereafter. 

Since the hydrogen perturber remains in its ground 
state Is during the collision, the wave function |?/>) of the 
system (atom+perturber), is the product of the wave func- 
tions \ip)(A) of the perturbed atom and that of hydrogen 



in its ground state |100)rm: 
IV) - M ( A)|100) {fl) 



(23) 



The eigenenergy of the hydrogen perturber in its ground 
state is denoted by Ejj. Following Anstee & O'Mara 
we use second-order time-independent perturba- 
tion theory for the eigenenergy of the system. 

We also use the unperturbed atomic wave functions 
defined by: 



H \nlmi) = E° lmi \nlmi); \ip) {A ) = \nlmi). 



(24) 



They form a complete set of eigenf unctions \nlmi) and 
eigenvalues E° lm[ . The perturbed eigenvalues of the sys- 
tem of Hamiltonian H can be written as: 



£ i = ^ + <i|V|i> + ^^ 



i\V\j)(J\V\i) 



(25) 



the summation being over all product states \j) = 
\^)(A)\^)(H) except \i) = \nlmi)( A )\ 100) (#), and where 



E 10 — TP 
a ' — a nlmi 



Quenching is neglected and thus we consider only the 
subspace nl (21+1 states) and we denote the product state 
of the two separated atom at R = oo by \Mi) (Anstee & 
O'Mara lTMTll . 



\Mi. 



| nlmi) w 1 100) 



(26) 



We introduce an effective potential V e t / , the diagonal el- 
ements of which arc defined as 



(M l \V e}f \Mi)=E Ml 



E Mi 



(27) 



The evaluation of the second order contribution to the 
above expression l|25|l is very difficult because it requires a 
summation over the complete set of eigenstates of the two 
isolated atoms. However, the expression can be greatly 
simplified by using the Unsold approximation, i.e. by re- 
placing E® — Ej by the constant value E p . With this ap- 
proximation, the expression for V e ff becomes 

{MAVefAMi) = (MMMi) - JL((A/i|v|M,)) 2 



1 

~E D 



E„ 



{Mi\V 2 \Mi), 



(28) 



where (Mi\V e ff\Mi) is the so-called Rayleigh-Schrodinger- 
Unsold (RSU) potential which is used throughout the 
ABO model. Then, for computing the RSU potential, ex- 
plicit expressions for the wave functions are needed. 

The perturbed atom is modelled as a single valence 
electron outside a spherical ionic core. The chosen ra- 
dial functions are of Coulomb type (Seaton I1958|l . The 
quantum defect and thus the effective quantum number 
n* = [2(Eoo — Eni)]^ 1 / 2 is a fundamental ingredient of 
the calculation of the interaction potential used in the 
ABO model. We refer to the ABO papers for more details 
concerning the interaction potential. 

The S'-matrix requires the time-dependent Schrodinger 
equation to be solved, i.e. 



{H + H P + V eff )U(t)) 



. d\m) 

1 dt 



(29) 



where the interaction potential V is now replaced by V e t / . 
The ^-scattering matrix is defined by: 



S = U (+oo, — oo) 



(30) 



where U(+oo, —oo) is the evolution operator in the inter- 
action representation, where 

${t + s) = U^(t + s,t)^(t)=e^ ffo+Hp >il;(t), (31) 

is the corresponding wave function. The wave function is 
expressed in terms of the basis formed by the eigenvectors 
\Mi) of (H a + Hp) given by (gBJ so that: 



|^(*)>=$>*«W e " 

Mi 



lE"t 



Mi); 



(32) 



$(£)) = Y^a Ml (t)\M l] 
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and we obtain a system of (21 ■ 
tial equations of first order: 



1) coupled linear differen- 



. da Ml (t) 
1 dt 



i 



= E a M k (t)(Ml\Veff\M k ) 



(33) 



M k = -l 

We have to solve this system for the initial conditions: 



0m, (-oo) =5(Mi,Mi), 

where S is the Kronecker symbol and 

(Mi\S\Mi) = {Mi\S\Mi) = a Ml (+oo). 



(34) 



(35) 



Thus the system has to be solved (21 + 1) times for obtain- 
ing all the elements of the S-matrix for the (nl) subspace. 

For the calculations, two frames of reference are used. 
The first one is the atomic frame in which the perturbed 
atom is stationary at the origin, the quantization axis is 
taken as perpendicular to the collision plane (b,v) and the 
second one is the rotating or molecular frame in which the 
interatomic axis is taken as the quantization axis. These 
reference frames are shown in Fig. This choice of frames 
gives the simplest relationship between states quantized 
relative to the atomic and the rotating frames respectively 
(RoueffQSZU. 

The interaction potential V has a cylindrical symmetry 
about the internuclear axis (OZ) and so V e ff is diagonal 
in the rotating frame. 

A transformation of the angular momentum states 
from the atomic basis to the rotating frame is achieved 
by using the rotations matrices fMessiah ll96Tfl . The Euler 
angles are given by (j3, f , f ) and equation l|3~3"| becomes 



. da Ml (t) 
1 dt 



E 

M k M! 



^Zm^mUv "M k (t) (M l | V eff | Mt) 



,t(0 



(36) 



when T> is the rotation operator. 

For a given quantum number I we obtain (21 + 1) cou- 
pled differential equations. These coupled equations must 
be solved for obtaining the (21 + l) 2 S-matrix elements. 

The interaction matrix is Hermitian and the S'-matrix 
is unitary and symmetric. If the quantization axis is per- 
pendicular to the collision plane the selection rules are: 
Am/ = ±2, ±4.... This rule will disappear when we take 
an average over all possible orientations of the collision 
plane, since the collisions are isotropic. 



6. Depolarization calculations for I = 1 

Let us consider collisions between an atom in an I = 1 
state and hydrogen in its ground state Is. The state of 
the system during the collision becomes explicitly: 

|V(*)> = ai(i)|ll)|00)e- i ^ t + a (t)|lO)|00)c- i£ o ' 

+ a-i(t)\l-l)\00)e- iE -^ (37) 



and the coupled differential equations become (Anstee & 
O'MaraEMJ: 



. da (t) 

1 m 

. da-i(t) 
1 dt 

where: 



a (t)(l\Veff\l) 

a_i(*)V+ +c.- 2il3 ai (t)V- 



(38) 



Vi = i(<l|K//|l)±(0|K//|0>). 



(39) 



The equations are integrated with a Runge-Kutta-Merson 
routine and then we obtain the transition matrix ele- 
ments in the \almi) basis for a given velocity and impact- 
parameter. 

We recall that we need, for depolarization calculations, 
the T-matrix elements in the \aJMj) basis. Since the spin 
is neglected, the S'-matrix is diagonal in S and its elements 
do not depend on Mg. Following equation (|20|) we obtain: 



(aJMj\T\aJM'j) 



E 



(_l)2S-M+M J+ M^ 2J + 1 ) 



mi,m' v Ms 

S I J 

M s mi -Mj 

S I J 



(40) 



1 M S m[ <«^,|T|aimJ) 

The matrix elements (aJMj\T\aJM'j) are obtained in a 
basis where the quantization axis is perpendicular to the 
collision plane. We have now to perform an average of the 
transitions probabilities over all possible orientations since 
collisions are isotropic. This angular average is denoted by 

( )av 

Using formula (16) by Sahal-Brechot (1973), it can be 
shown that, in the irreducible tensorial operator basis, the 
angular average of the depolarization transition probabil- 
ity is given by: 

(P k (aJ,b lV )) av = J2\( aJ v\T\a J m')| 2 

- ^2 ( a J M|7> J M') (" J v\T\a J v')* (41) 



E(-d 2J 



+k+n— fi' 



J J k \ f J J k 



-v' [»! X I \v -fx -x 



Owing to the selection rules on the 3j-coefficients, the 
summation over \ 1S reduced to a single term, since x — 

-( m '-^) = -(m-^)- 

The depolarization cross section may then be com- 
puted from: 



a k (aJ,v) = 2tt / (P k (aJ,b,v)) av b db 



(42) 



At small impact-parameters, the (P k (aJ,b,v)) av values 
become highly oscillatory, and it is impractical to integrate 
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numerically. In addition, at these separations the pertur- 
bation theory breaks down and the calculated RSU poten- 
tials become invalid. Thus, an impact-parameter cutoff is 
used (Anstee & O'Mara lTMTjl : 



r OO 

<r k (aJ, v) ~ irbl + 2tt (P k (aJ, b, v)) av b 

J b 



db 



(43) 



Integration over the velocity distribution for a tempera- 
ture T can be performed to obtain the depolarization rate 
which is given by: 



D k {aJ,T) ~ n H J vf(v)dv\wb% 

(P k (aJ,b,v)) av bdb 



(44) 



+2tt 



bo 



3a 



where bo is the cutoff impact-parameter and here bo 
as in Anstee & O'Mara lfTMT|) (see Sect.©. 

7. Results 

7.1. Our results for 1 Pi and 2 Pz/2 levels 

The perturbed atom is modelled as a single valence elec- 
tron outside a spherical ionic core. If the energy of the 
state \al) of the valence electron is denoted by E a i and 
the binding energy of the ground state by E^, the bind- 
ing energy of the level \al) is (E^ — E a i) and is re- 
lated to the effective principal quantum number by n* = 

[2(£oo-£n0r 1/2 - 

For example: consider the case of the neutral cal- 
cium line 4227 A (transition 4s 2 — 4s4p). The energy of 
the 4p level and the binding energy are ~ 0. 1077a. u.(^ 
23652 cm" 1 ) and ~ 0.2245a.u.(49305.72 cm" 1 ) respec- 
tively (Wiese et al I1969[) . Hence the effective principal 
quantum number is n* ~ 2.077. 

Note that if the parent configuration of the level is not 
the same as that of the ground state, this should be ac- 
counted for by using the appropriate series limit for the 
excited parent configuration instead of the ground state 
binding energy. A particular example was explained in 
detail by Barklem, Anstee & O'Mara (1 998(1 (see also 
Barklem, Piskunov & O'Mara I2000|) . As in Anstee & 
O'Mara (1991), we can calculate the cross sections as a 
function of n* . Table gives a k (al,v) and a k (a^,v) as a 
function of n* for a typical relative velocity v — 10 km s . 

We thus can obtain the depolarization cross section for 
any level of any atom by interpolation of the values given 
in Table for the appropriate n* value for the state of 
interest. We notice that the dependence of cross sections 
on velocity is very close to 



a k {aJ,v)=G k {aJ,vo){-)- xk{aJ \ 



(45) 



where vo is the velocity at which the cross section is known 
(10 km s _1 ). Table|31shows the so-called velocity exponent 
X k (aJ) (Anstee & O'Mara HWfy . In order to obtain the 



n* 


a L {al) 


<7 2 (al) 


a>§) 




<x a (af) 


1.5 


203 


192 


79 


141 


134 


1.6 


233 


228 


90 


165 


157 


1.7 


276 


258 


100 


187 


176 


1.8 


323 


301 


114 


217 


203 


1.9 


377 


349 


130 


252 


233 


2 


438 


405 


147 


290 


268 


2.1 


507 


466 


167 


334 


307 


2.2 


581 


533 


189 


381 


349 


2.3 


663 


605 


212 


432 


393 


2.4 


751 


683 


237 


487 


442 


2.5 


849 


769 


266 


549 


496 


2.6 


949 


869 


295 


616 


562 


2.7 


1058 


970 


327 


685 


627 


2.8 


1180 


1075 


362 


761 


691 


2.9 


1299 


1176 


396 


835 


753 


3 


1423 


1282 


431 


911 


817 



Table 2. Variation of the depolarization cross section, for 
a relative velocity of 10 km s _1 , with effective principal 
quantum number. Cross sections are in atomic units. 



n* 


X l (al) 


A z (al) 


V(a|) 


A>f) 


A>§) 


1.5 


0.269 


0.256 


0.201 


0.246 


0.240 


1.6 


0.270 


0.257 


0.207 


0.246 


0.234 


1.7 


0.265 


0.247 


0.211 


0.244 


0.235 


1.8 


0.259 


0.249 


0.213 


0.242 


0.240 


1.9 


0.256 


0.242 


0.215 


0.240 


0.234 


2 


0.251 


0.235 


0.216 


0.237 


0.235 


2.1 


0.253 


0.237 


0.220 


0.235 


0.218 


2.2 


0.248 


0.235 


0.219 


0.232 


0.223 


2.3 


0.246 


0.235 


0.221 


0.231 


0.221 


2.4 


0.242 


0.231 


0.220 


0.230 


0.225 


2.5 


0.240 


0.224 


0.221 


0.228 


0.224 


2.6 


0.242 


0.227 


0.223 


0.228 


0.220 


2.7 


0.241 


0.223 


0.224 


0.228 


0.226 


2.8 


0.244 


0.226 


0.228 


0.230 


0.223 


2.9 


0.249 


0.223 


0.221 


0.230 


0.208 


3 


0.243 


0.229 


0.229 


0.2208 


0.225 



Table 3. Velocity exponents X k (aJ) corresponding to the 
cross sections of Table 2. 



depolarization cross section for any velocity v from Tables 
|2| and [3] equation l|45l) can be used. 

Figures and 0] display the depolarization cross sec- 
tions (computed directly) and the corresponding depo- 
larization rates for the Cal Ap 1 Pi{n* ~ 2.077), Mgl 
3p 1 P 1 (n* ~ 2.030) and the Nal 3p 2 P 3 / 2 (n* ~ 2.117) 
states. 

To obtain the depolarization rates we must integrate 
over the Maxwell distribution of velocities. Using equation 

(I44|l , it can be shown that for a mean velocity v = \J^f^, 

the depolarization rates can be expressed as (Anstee & 
O'Mara UMB 
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Fig. 3. Depolarization cross sections as a function of rel- Fi §- 4 " Depolarization rates as a function of temperature, 
ative velocity. 



D k (aJ,T) = (-)^ xk ^r(2-^X k (aJ))v a k (aJ,v ) 

IT Z 

x (l)l-A fc (aJ) (46) 



Since f(v) depends on the reduced mass /i, the depolar- 
ization rate is specific to a given atom. In order to obtain 
depolarization rates for any level of any atom by simple 
interpolation we applied the approximation: /i = run, tuh 
being the hydrogen mass. We verified that this approxima- 



10 



M. Derouich et al: Semi-classical theory of collisional depolarization 



n* 


T=5000K 


T=6000K 


D'(al)/n H 


D*(a§)/n H 


D\aY)Jn H 




1.5 


0.5338 


0.3931 


0.5705 


0.4206 


1.6 


0.6135 


0.4477 


0.6556 


0.4789 


1.7 


0.7188 


0.5220 


0.7679 


0.5585 


1.8 


0.8363 


0.6047 


0.8944 


0.6470 


1.9 


0.9716 


0.7000 


1.0395 


0.7492 


2. 


1.1260 


0.8084 


1.2050 


0.8653 


2.1 


1.2971 


0.9289 


1.3885 


0.9944 


2.2 


1.4844 


1.0599 


1.5891 


1.1349 


2.3 


1.6875 


1.2021 


1.8066 


1.2872 


2.4 


1.9043 


1.3563 


2.0390 


1.4523 


2.5 


2.1452 


1.5270 


2.2982 


1.6352 


2.6 


2.4065 


1.7086 


2.5774 


1.8297 


2.7 


2.6814 


1.9010 


2.8709 


2.0348 


2.8 


2.9634 


2.1015 


3.1695 


2.2469 


2.9 


3.2665 


2.3097 


3.5001 


2.4707 


3. 


3.5842 


2.5237 


3.8410 


2.6992 



Table 4. Depolarization rates for I = 
S = and J = 1; S = § and J = § 
T=6000K, as a function of n*. Depolarization rates are 



1. Each column: 
for T=5000K and 



given in 10 14 rad. m 3 s 



tion introduces a less than 4% error in the depolarization 
rates. 

Table 0| gives the alignment depolarization rates as a 
function of n* for the temperatures of 5000 K and 6000 
K. 

7.2. Comparison with other methods 

In Figure[S]we compare our alignment depolarization rates 

with quantum chemistry depolarization rates (Kerkeni 

2002 ) and the alignment depolarization rates obtained by 

replacing the RSU potential by the Van der Waals poten- 
cy 

tial, V — where Cq is the Van der Waals constant. 

We display only the k = 2 case which is related to 
the linear depolarization (alignment). Our results show 
rather good agreement with quantum chemistry calcula- 
tions. However, the Van der Waals potential underesti- 
mates the depolarization cross section. To explain this we 
now examine the sensitivity of the depolarization cross 
section to the RSU potential. 
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Fig. 5. Depolarization rates for k — 2 as a function of 
temperature. Full lines: our results; dotted lines: quantum 
chemistry calculations (Kerkeni 2002); dot-dashed lines 
Van der Waals approximation. 
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8. Dependence of depolarization cross sections on 
potentials 

As a first check on the sensitivity of the results to the na- 
ture of the potential we introduce a localised perturbation 
to the interaction potential. We have chosen the case Cal 
4p for the check. The interaction potential is multiplied 
by a Gaussian magnification factor of the form (Anstee & 
O'Mara nMT]) 



G{R) = l + exp{-{R- R ) 2 ). 



(47) 



The (7-symmetry interaction (ls,4pa) is more impor- 
tant than the 7r-symmetry interaction (Is, 4p7r). Therefore 
only the (ls,4pcr) interaction has been altered. Rq is the 
position of the maximum perturbation. 

Fig. shows the effect of the Ro variation on the de- 
polarization cross section for three different velocities. We 
remark that the enhancement effect on the cross section 
shows a peak in the intermediate range (10ao < Rq < 
18ao). We conclude that the depolarization cross section 
depends most strongly on the intermediate region of the 
potential. It is for this reason that the Van der Waals inter- 
action potential, which is inaccurate in this intermediate 
region, underestimates the depolarization cross sections. 

To study the sensitivity of the cross section to close 
collisions, a second check is made. We have calculated the 
alignment cross sections by varying the cutoff &o in equa- 
tion l|4H|l. Fig. [3 shows the effect of the cutoff variation at 
v — 10 km s _1 and v = 30 km s _1 . For bo < 8ao the cut- 
off variation does not have any effect on the cross section 
calculations. 

Therefore it can be concluded that short-range col- 
lisions which involve ionic and exchange interactions do 
not influence the depolarization cross sections and hence 
the depolarization rates. The principal differences between 
the RSU potentials and those from quantum chemistry oc- 
cur at small interatomic separations. Since the important 
impact-parameters are at medium interatomic distances, 
this explains why we obtain a rather good agreement be- 
tween the ABO and quantum chemistry results. However, 
a more detailed study would be desirable but is outside 
the scope of the present paper because the exact quantum 
chemistry potentials are needed for quantitative compar- 
ison. 
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Fig. 6. Depolarization cross section enhancement due to 
a Gaussian local perturbation of the potential. Cross sec- 
tions are calculated for the relative velocities: 3 km s -1 , 
10 km s _1 and 17 km s _1 . 
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Fig. 7. Behaviour of the depolarization cross sections from 
equation (|43|l with the cut-off &o- For a relative velocities 
v = 10 km s _1 and v = 30 km s _1 . 



9. Conclusion 

In conclusion, we have shown that the ABO theory is 
a powerful tool for quickly calculating a great number 
of depolarization cross sections for collisions with atomic 
hydrogen in the physical conditions of the solar atmo- 
sphere. The present paper is limited to p states for testing 
the method. Comparison with quantum chemistry results 
gives us the conviction that it would be useful to extend 
our calculations to higher ^-values which are also of im- 
portance for solar line polarization. Although less accu- 
rate than quantum chemistry calculations, an extension 
to heavy atoms and higher Z-values is easier and faster 



within the present semi-classical method. We will also ex- 
tend our semi-classical theory to depolarization cross sec- 
tions by collisions with singly ionized atoms. This will be 
the subject of further papers. 

Acknowledgements. This paper is dedicated to Jim O'Mara's 
memory, who followed with interest the progress of the present 
work and suddenly passed away on April 27, 2002. The authors 
thank the referee for valuable suggestions which improved the 
presentation of this work. We are indebted to B. Kerkeni who 
has communicated a number of her results before publication. 



12 



M. Derouich et al: Semi-classical theory of collisional depolarization 



References 

Ackermann H., 1966, Colloques Internationaux du CNRS 

(Paris, France: editions du CNRS), N° 164, 82 
Anstee S.D., & O'Mara B.J., 1991, MNRAS, 253, 549 
Anstee S.D. PhD thesis, Univ. Queensland, 1992 
Anstee S.D., & O'Mara B.J., 1995, MNRAS, 276, 859 
Anstee S.D., O'Mara B.J., & Ross J.E., 1997, MNRAS, 284, 
202 

Barklem P.S., & O'Mara B.J., 1997, MNRAS, 290, 102 
Barklem P.S., & O'Mara B.J., & Ross J.E., 1998, MNRAS, 
296, 1057 

Barklem P.S. PhD thesis , Univ. Queensland, 1998 
Barklem P.S, Anstee S.D., O'Mara B.J., 1998, PASA 15(3), 
336. 

Barklem P.S., Piskunov N., & O'Mara B.J., 2000, A&A, 142, 
467 

Bommier V., & Sahal-Brechot S., 1982, Solar Physics, 78, 157 
Bommicr V., & Molodij G., 2002, A&A, 381, 241 
Brueckner K.A., 1971, ApJ, 169, 621 

Fano U., & Racah G., 1959, Irreducible Tensorial Sets (New 

York; Academic Press) 
Fano U., 1963, Phys. Rev., 131, 259 

Gangrsky Yu.P., Karaivanov D.V., Marinova K.P., Markov 
B.N., Melnikova L.M., Mishinsky G.V., Zemlyanoi S.G., 
& Zhemenik V.I., 1998, Eur. Phys. J. A, 3, 313 

Gaupp A., Kuske P., & Andra H.J., 1982, Phys. Rev., A26, 
3351 

Gordeyev E.P., Nikitin E.E., & Ovchinnikova M. Ya., 1969, 

Can. J. Phys., 47, 1819 
Kerkeni B., Spielfiedel A., & Feautrier N., 2000, A&A, 358, 373 
Kerkeni B., 2002, A&A, 390, 791 
Kerkeni B., & Bommier V., 2002, A&A, 394, 707 
Messiah A., 1961, Mecanique Quantique (Paris: Dunod) 
O'Mara B.J., 1976, MNRAS, 177, 551 
Omont A., 1977, Prog. Quan. Elec, 5, 69 
Roueff E., 1974, J.Phys. B, 7, 185 
Roueff E., & Suzor A., 1974, J.Phys. B, 35, 727 
Sahal-Brechot S., 1973, A&A, 32, 147 
Sahal-Brechot S., 1977, ApJ., 213, 887 
Seaton M.J., 1958, MNRAS, 118, 504 
Stenfio J.O., & Keller C.U., 1997, A&A, 321, 927 
Unsold A.L., 1927, Zeitschrit fur Physik, 43, 574 
Unsold A.L., 1955, Physik der Stern Atmospharen (Zweite 

Aufiage) 

Wiese W. L., Smith M. W., & Miles B. M., 1969, Atomic 
Transition Probabilities - Sodium through Calcium. U. 
S. Dept. of Commerce, National Bureau of Standards, 
Washington 



